Universality of the conductance in quantum dot transport 
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Abstract 

We revisited the scaling behavior of the transport properties of a quantum dot system described 
by the spin-1/2 Anderson model using analytical methods. In the low temperature limit we show 

<D . 

that the conductance has a universal behavior with universality between temperature and bias. We 
compare this result with the empirical formula used to fit the experimental data for conductance in 
the case of the equilibrium transport through a single channel quantum dot. In the high temper- 
ature limit the conductance obtained from the Anderson model is compared with previous results 

o : 

obtained from the Kondo model. The universal behavior is present also in the high temperature 
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(N 
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o 
o> 
o 
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limit. These results are in good agreement with the Renormalization group calculations. 
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I. INTRODUCTION 



Recent experimental results demonstrated the importance of the electronic correlations 
for the behavior of the low temperature transport of quantum dots and break junctions 
devices containing C 60 . A quantum dot (QD) consist of a confined region of electrons of 
size I ~ Xb, (Ab - the de Broglie wave-length), coupled to leads via tunnel barriers. The 
conduction of such a device was described by the well known Anderson model l|, |2j, in terms 
of the interacting constants A = 7riV(0)| Vkd\ 2 ( Vkd - the interaction between the conduction 
electrons and ci-impurities) and Coulomb interaction U. Most of the transport studies for 
such a configuration are carried in the T = limit, in the case of strong (U/A ^> 1) and weak 
(U/A <C 1) coupling between the conduction and localized electrons. In the strong coupling 
and zero temperature limit the Coulomb interaction U provides the largest energy scale. 
Although the T = case provides a reasonably answer for the transport problem, the case 
of a finite temperature (T ^ 0) should be carefully considered as the transport properties 
of the system can be significantly influenced by thermal effects. For A < T < U quantum 
fluctuations are small and transport is dominated by charge effects. This regime is well 
described by the concept known as Coulomb blockade and in this situation the temperature 
dependence of the conductance G(T) was calculated by Beenaker [3j. At low temperature the 
conductance presents a series of equidistant peaks as function of gate voltage V separated 
by a spacing U. These peaks correspond to a fractional number of electrons in the dot 
and separate the Coulomb blockade valleys corresponding to either an even or odd number 
of electrons. In the regime T « A « U the strong quantum fluctuations can lead to 
important modifications of the Coulomb blockade scenario. For an odd number of electrons 
in the quantum dot (called now a Kondo quantum dot) the system can have a net spin 1/2 
and the Kondo effect can develop. As a consequence the conductance in the odd electron 
valleys is enhanced, turning those valleys into plateaus of near perfect transmission. 

The experimental measurement of the system's conductance at temperatures below the 
Kondo temperature (Tk) have been a difficult task for physicists. In particular, the low tem- 
perature regime defined above (T <C A) requires a direct control over the interaction between 
localized and conduction electrons. Recent advances in nanotechnologies allow a better con- 
trol over the geometry of the system's electrodes and implicitly over the electron interactions. 
The temperature dependence of the system's conductance 0, 5, (| have been fitted using the 
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empirical formula G EK (T) = G [l + 2 1 ' S - 1 {T /T K ) 2 ]- S where T K = T K /{2 l ' 2s - l) 1 / 2 with 
s = 0.22 and Gek{Tr) = Go/2. On the other hand, in the low temperature limit, theo- 
retical estimations of the system's conductance lead to G(T) = Go[l — tt 2 (T/Tk) 2 /16] with 
Go — 2e 2 /h; this result is given by the temperature dependence of the Kondo resonant peak 
in the system's density of states. Although there is a good agreement between experimental 
and theoretical results in the very low temperature limit (T <C Tk), in the limit of higher 
temperatures experimental and theoretical data behave differently. 

Here, we revisit the problem of universal behavior of the transport properties in quantum 
dot systems. In the equilibrium case the scaling behavior was obtained using the microscopic 
Anderson model with a finite on-site Coulomb interaction U in the non-crossing approxi- 
mation [7] , but this method can give only qualitative results. On the other hand, the origin 
of the differences between the theoretical conductance G(T) obtained in the Fermi-liquid 
approximation and the experimental conductance Gek was not discussed in the literature. 
In Section II we present the calculation of the system's conductance in the equilibrium case 
for finite on-site Coulomb interaction U and discuss the possible origin of the enhancement 
in the Kondo temperature. The universal behavior for the non-equilibrium case is presented 
in Section III using the microscopic Anderson model. This problem was treated first by 
Appelbaum [§| and Anderson ^ and more recent reconsidered by Schiller and Hersfield 10] 
using the Touluse exact model. The same situation, but for the case of low and high bias 
voltages, was reconsidered by Oguri using the Yamada [/-expansion [141 ] . Using a 

similar model we reconsider the problem in connection with recent experimental results ob- 
tained by Grobis et al. 21( and recalculate the scaling function of the system's conductance. 
In Section IV we discuss our results in connection with other results from the literature. 



II. UNIVERSAL SCALING IN EQUILIBRIUM TRANSPORT THROUGH A 
QUANTUM DOT 

In this section we will present a detailed analysis of the universal behavior of a single 
channel quantum dot using the Anderson model. We will consider both low temperature 
(T <C Tk) and high temperature (T ^> Tk) regimes. We will discuss the possible origin of 
the differences between the theoretical and empirical fits to the system's conductance in the 
low temperature regime. 
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A. Model and Green functions 



We start from the simple symmetric Anderson model described by the Hamiltonian: 

H = H iead + H dot + H int (1) 



where: 



Hlead — tkc\o-Cka , (2) 



H dot = e «*4 d <r + U (d\ d i ~ ^ \4 d i- ~ ^) ~ J 



(3) 



and 



Hm = + M • ( 4 ) 

Act 

This Hamiltonian has been used by Anderson to describe the occurrence of the magnetic 
moments in dilute alloys. The first part of the Hamiltonian, Hi ea d, describes the conduction 
electrons in the leads in terms of the creation (c k(7 ) and annihilation (c^) operators. The 
second part of the Hamiltonian, Hdot, similar to the Anderson's model d-localized level, 
describes the electrons localized in the quantum dot in terms of the operators S a and d a . 
Finally, H in t describes the interaction between the conduction electrons and the localized 
electrons. 

In order to calculate the physical quantities described by this Hamiltonian we will use 
the Green function: 

G^^iun) = GQ X (iu n ) - S(zw n ) , (5) 

where uo n = (2n + l)irT. Go(ioo n ), the bare Green's function in the absence of the Coulomb 
on-site interaction U, can be calculated as: 

Safari) = ' ( 6 ) 

iu n - e d + i/\sgnuj n 

where A = irN(0)\Vkd\ 2 with N(0) being the density of states at the Fermi level for the 
conduction electrons. The self-energy S(io; n ) describes the effect of the on-site Coulomb 
interaction and can be evaluated using the Yamada perturbative [/-expansion [14j to the 
(u 2 , T 2 ) order in the region U/(ttA) ~ 1 and < e^/A < 1. The energy scale of the model, 
for the U = case, is given by the Kondo temperature Tk'- 

^ = Aexp(^), (7) 



where E d = e d + Aln(D/A)/7r, and D is the bandwidth. Note that Tk given by this 
expression has to be considered only as an energy scale, a better expression for the Kondo 



temperature, which can be used to fit the experimental result, being T r K = y/ ATk- 

The density of states corresponding to the localized electrons in the quantum dot can be 
obtained from their Green function as: 

' W ' T ) = ~Z /mG( ^' T) = Z 7 7-1 / -J',*, ^2 , (8) 

Ti 71 [uj - E d (uj,T)\ +[/mS(w,r)] 

where E d (uo, T) = + ife£(a>, T). In the following we will consider two different situations, 
the case of low temperatures, (uj,T) « Tk, and the case of high temperatures, (00, T) » 



1. Low temperature regime 

In this case ReYi(uj,T) ~ — e d and /mS(w,T) = A + Ei(u;,T) where £i(a>) is given by: 



Si( W) T) = aA^^L (9) 



cu 2 + (ttT) 2 

The value of the eo„ 8tant „ = , V 16 was given by Costi e, aL M , for the s ym ,„e t nc 
Anderson model using the Yamada [14J] expansion of the se 
value obtained in the Fermi liquid approximation by Nozieres 
constat a, the resulting temperature dependence of the system's conductance describes the 
experimental results only qualitatively. An analytical expression of the constant a from the 
self-energy, which can change quantitatively the Kondo temperature, was given by Kirchner 



f-energy and is exactly the 
Ji|. Using this value for the 



et al. 



17|: 



7T 4 I^ 2 R-l \e d \ 



(10) 



8e (3/2+C7) sin^TTn^) D ' 

where C = 0.57721 is the Euler constant, W = 0.4128 is the Wilson number, R = 2 
is the Wilson ratio, and 2D is the bandwidth of the conduction electrons. The important 
parameter of the model is the product aA which for the Fermi liquid model was calculated as 
aA = 0.0239 D. Using the result given in Ref. [13] , taking for \e d \/D = 0.81 and A/D = 0.2 
we obtained aA = 0.0236, which is a reasonable value compared with the value from the 
Fermi liquid model. We can get information about the magnitude of this parameter from 
the experimental data obtained for the temperature dependence of the Kondo resonance 
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width [18(1 . Based on Eq. ([9} the imaginary part of the self-energy can be written as: 



ImE(T,u) = A 



1 + 



T, 



K 



+ 



ttTV 



(11) 



where T£ = T K /y/a. Using the above expression the electronic density of states correspond- 
ing to the QD electrons can be approximated as: 

-l 



p(T,u) 



1 K 



1 K 



(12) 



or if we introduces the notation 7 2 (T, T K ) = (vrT) 2 + (T K ) 2 , as: 



p(T,u) 



7T 



A 



1 + 



ttTV 



1 2 (T,T* K ) 
7 2 (T,T^)+^ 



(13) 



This form for the density of states is specific for the Fermi liquid behavior and describes the 
low temperature domain defined as (u),T) « Tk- 



2. High temperature regime 



In the high temperature regime, u 3> Tk, the effects of the Kondo state on the system's 
conductance can be evaluated using the electronic self-energy in terms of the t-matrix ap- 
proximation, £(co>) ~ t(u). The electronic t-matrix can be calculated using the Anderson 
model along with an appropriate decoupling in the equation of motion, a method which is 
equivalent to the calculation of the t-matrix with respect to a perturbation expansion up 
to the sixth order in the interaction term V^d 19|- The connection between the electronic 



t-matrix and the <i-electrons density of states was discussed first by Langreth [20] and is 
given by 



Im t(u) 



N(0) 



p{u,T) . 



(14) 



In the high temperature limit, for a correct description of the system's properties, the 
t-matrix approximation has to include effects related to the inelastic scattering of the con- 



duction electrons on t 
Mamada and Shibata 
can be done in the U - 



re spins attached to the G?-electrons in the quantum dot. Following 



19] the evaluation of the t-matrix for the symmetric Anderson model 



oo limit: 

Im t(u) = -7rJV(0) 



(15) 
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where the running coupling constant J(uS) = 1/ ^7riV(0) In in the high temperature 
(\u\ » T K ) regime. Based on the above equations the electronic density of states can be 
approximated as: 

p(u,T) ~ 4" T7: ■ ( 16 ) 

A 16 In Tjl- K J 

B. Conductance 

The conductance is one of the most investigated physical quantities in quantum dot 
systems. For the transport properties of a single quantum dot system we expect a universal 
behavior of the conductance G(T/Tk), behavior which is specific to the Kondo problem. 
However, the experimental facts prove that there are deviations from the specific Fermi 
liquid behavior in the low temperature regime, T <C Tk- As a consequence the experimental 
data were fitted with an empirical formula Gek(T). In the following we will calculate the 
conductance behavior for this regime using the Anderson model. On the other hand the high 
temperature behavior is obtained taking the results from the perturbation theory applied to 
the Anderson model. 



1. Low temperature regime 

To obtain the temperature dependence of the system's conductivity in the low tempera- 
ture regime we start with the general equation for the current: 

I(V,T) = — J dup(u)[f(u-eV)-f{u + eV)], (17) 

where V is the applied external bias and f(uo) is the Fermi function. The integral in the 
above definition can be calculated and using the definition of the system's conductance, 
G(V,T) = dI(V,T)/dV, we find: 

where 7 = 7(T, T£) and ^I^(z) = d^?(z)/dz is the trigamma function. Using the asymp- 
totic behavior of the trigamma function, ^f^(z) ~ l/z + l/2z 2 , in the low voltage limit 
eV/k B T K <^ 1, the system's conductance can be approximated as 

G(0,T)~G 



(18) 



2 (5)' 



(19) 
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a result which is in good agreement with Fermi-liquid behavior, i.e., the conductance de- 
creases quadratically with temperature T in the low temperature regime, T « T^. This 
expression for the system's conductance is very similar to the empirical formula used in 
fitting the experimental data Gek(Q,T) = Gq [l — ct (T/T^) 2 ], where ct is determined by 
the definition of the Kondo temperature Gek(T k ) = G /2. 



2. High temperature regime 



In the high temperature regime, o> >> T K , the system's conductance G(T) can be eval- 

221: 



uated using the general relation 

G(T) = Go 



du 



df(T) 
du 



[-N(0)lmt(u)] 



(20) 



a relation which combined with the expression for the electronic t-matrix leads to: 

3tt 2 1 



G(T) ~ G - 



16 In 2 ^ ' 



(21) 



where the J 3 (u;)-term has been neglected. This result was obtained first by Kaminski et al. 



24] from the Kondo Hamiltonian. We mention that in the recent analysis by Borda et al. 
23 ] the inelastic scattering gives a 1/\b. 2 (uj/Tk) term with the increasing of the parameter 



U/A. This conclusion is in agreement with the above calculation of the conductance since 
the Kondo model is equivalent to the Anderson model in the limit of large U / A and u « U. 



III. UNIVERSAL SCALING IN NON-EQUILIBRIUM TRANSPORT THROUGH 
A QUANTUM DOT 

In this section we analyze the nonequilibrium transport properties of a quantum dot 
system. We will consider in detail the low temperature regime, T <C T%, using the same 
formalism as in the equilibrium case to identify the possible universal behavior of the sys- 
tem's conductance, occupancy of the quantum dot, and magnetic spin susceptibility in the 
nonequilibrium case. 



S 



A. Conductance 



Following the same procedure as in the equilibrium case, we will need to evaluate first 
the electronic density of states for our system. The electron Green's function for the 



nonequilibrium case was ca. 



Renormalization theory llj by Oguri as 



culated in the e 



G r (u) 



ectron-hole symmetric limit using the Hewson- 



(22) 



where A = zA. The parameters in the above equation are given by z = 1 — (3 — tt 2 /A)u 2 , 



u = U/(irA), and U = z 2 T ^(0,0), T ajCT i(ui, uj) being the renormalized vertex function. 
Accordingly, the system's density of states defined as p(ui,T) = —ImG r (uj)/Tc can be ap- 
proximated as: 



1 



ttA 1+ i (i)V 



(23) 



±(eV) 2 + (nTy] 

Following the same procedure as in the equilibrium case, in order to evaluate the system's 
conductance we will start from the estimation of the nonequilibrium current in the presence 
of an external voltage. Starting from the general definition we find 

tanh(^) -tanh(-"- eV ^ 



I(V) 



duo- 



2T 



1 + a[u 2 + \{eV) 2 + (wT) 2 ] 



where a = 1/(2T )(U /irT) 2 . This integral can be performed exactly and we obtain : 



I(V) = —Imm 
aa 



1 a + ieV 

2 + 2%T 



(24) 



(25) 



where *ff(x) is the digamma function and a 2 = 1 + - + §(eV^) 2 + (7rT) 2 . The imaginary part 
of the di-gamma function can be evaluated and the current I(V) becomes: 
2e V2T K 



I(V) 



h 



1 + 8 ( T K ) + 2 \T K 



arctan ■ 



eV 
V2T R 



3 I eV_ \ I 7T 2 It 
T K ^ 2 \T K 



1/2 



(26) 



or in the low voltage limit (V — > 0) 

I{V) 



h 



l 



8^T K ) ' 2\T K ' 

Accordingly, the system's conductance exhibits a scaling behavior of the form [G(T, 0) — 
G(T, V)]/G = F(eV/T K ,T/T K ) with 



; t^'t^)~8 \r K ) " i6 vtJ v?v 



(28) 



a result which is qualitative agreement with Ref. 



B. Particles number and magnetic susceptibility 

Additional parameters show an universal behavior in quantum dot systems. For example, 

for the single dot system the occupancy (number of particles) of the main quantum dot and 

the magnetic susceptibility of the system can be shown to present an universal behavior. 

Let consider for the beginning the occupancy of the main dot of the system. By definition 

2 r°° 

n d = / duf(u)ImG d (uj) , (29) 

where f(u) represents the Fermi function and G r d (u) is the electronic Green's function. The 
integral on the right hand side of the above equation can be performed and we get: 

n d = = , (30) 

where we used the Wilson ratio R defined as R = 1 + U/(nA). The system's magnetic 
susceptibility x(T) can be calculated as: 

X(T) = jT du,p{u, - /(«)] • (31) 

The calculation is relatively simple and we obtain 

X(T, V) = ,„,„ r ...... ■ P2) 



(f)+- 2 (i) 

Note that for A = Tk and R = 2 we obtain the scaling behavior in the Kondo regime 



IV. SUMMARY AND CONCLUSIONS 



We studied the universal behavior of a single quantum dot system at low temperatures 
described by the microscopic Anderson model using the many body methods. The main 
results can be summarized as follows: 

• the temperature dependence of the system's conductance in the equilibrium case is 
given by the second order corrections in the Coulomb interaction by an effective 
temperature dependent density of states. In the low temperature approximation, 

10 



T « Tk, the correction has the form G(T/T^) ~ a(T/T£) 2 , where the constant 
a depends on the model's parameters. This result can be a reasonable explanation 
for the existence of quantitative differences between the experimental data and the 
formula obtained using the standard Fermi liquid. The empirical formula Gke{T) 
does not contains a new physical picture, but new numerical coefficients which may 
be given by a more realistic model. 

in the high temperatures limit, T » T%, the system's transport properties in the 
equilibrium case were obtained using the t-matrix approximation for the Anderson 
model. The system's conductance was calculated as G{T/T* K ) ~ 1/ \n 2 (T/Tt), a re- 



sult obtained previously in the literature using the Kondo Hamiltonian |24j. This 
result shows that in the higher order of the perturbation theory we can obtain, us- 
ing a microscopic Hamiltonian, the same result as from the phenomenological Kondo 
Hamiltonian. 

the scaling behavior of the system's conductance in the equilibrium case was analyzed 
in literature using the Hewson-Renormalization procedure for the Anderson model. In 
our work, we obtained results similar to the proposed phenomenological description 
of the universal scaling in non-equilibrium transport (Grobis et al. [2l{|) using the 
microscopic model and many-body methods. These results reduces to the equilibrium 
case for the zero bias case. 
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